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ABSTRACT 

We analyze the physical properties of stellar clusters that are detected in 
massive star-forming regions in the MYStlX project—a comparative, multiwave¬ 
length study of young stellar clusters within 3.6 kpc that contain at least one 
0-type star. Tabulated properties of subclusters in these regions include physical 
sizes and shapes, intrinsic numbers of stars, absorptions by the molecular clouds, 
and median subcluster ages. Physical signs of dynamical evolution are present 
in the relations of these properties, including statistically significant correlations 
between subcluster size, central density, and age, which are likely the result of 
cluster expansion after gas removal. We argue that many of the subclusters iden¬ 
tified in Paper I are gravitationally bound because their radii are significantly less 
than what would be expected from freely expanding clumps of stars with a typ¬ 
ical initial stellar velocity dispersion of ~3 km s“^ for star-forming regions. We 
explore a model for cluster formation in which structurally simpler clusters are 
built up hierarchically through the mergers of subclusters—subcluster mergers 
are indicated by an inverse relation between the numbers of stars in a subcluster 
and their central densities (also seen as a density vs. radius relation that is less 
steep than would be expected from pure expansion). We discuss implications of 
these effects for the dynamical relaxation of young stellar clusters. 
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Introduction 


Stars form from hierarchically collapsing molecular clouds, which leads to clustered 
star-formation often reflecting the structure of the molecular cloud. Observations of massive 
star-forming regions (MSFRs) at different stages of their star-forming lifetime (typically <5- 
10 Myr), reveal diverse stellar cluster structure, which provides information about cluster 
formation, cluster dynamics, and future cluster survival or destruction— problems that have 


wider implications for the o rigin of stellar populations in galaxies (e.g., iLada fc Ladal 12003 


Portegies Zwart et al.ll2010l) . Important questions on these topics include: 


How are stellar cl usters formed? Do they freeze out of molecu lar clouds in a single 


crossing time fiElme green! l2000l: I Vazquez-Semadeni et ah! 12003 ) ? Ar e they built up 


by gradual star form ation over man y crossing times flTan et ah 
from the inside out fjPfalznei 


200611 ? Do they form 


201 ill? Or, a r e clusters formed from the hierarchical 


merger of subclu sters flMcMillan et ah! l2007l: iMaschberger et af 

20081: lBatell200^ ? 


2010: Bonnell et ah 


What is the origin of stellar clusters with different pr operties? Do stellar clusters 
start from extremely dense states with r\j 10® stars pc ^ flBatell2009al: iMoeckel fc Bate! 


2010l) ? Are observed differences in Gala ctic open cluster properties due to cluster evo¬ 


lution fjPfalznerll2009l: iGieles et al.ll2012[l ? Is the mass segregation that is seen in some 
regions primordial, indicating different star-for mation mechanisms for high and low 


mass stars, or can it be explained dynamically (IBonnell fc Davieslll998l: lAllison et ah 

2OO9I1 ? 


What environments lead to bound stellar clusters after the end of star formation? 

What roles do star-formation efficiency, cluster dynami cs, star-formatio n feedback ef¬ 
fects, and the tidal disruption of molecular clouds play fjKruiissenl l2012l) ? 


Advancements in multiwavelength st udies of young stellar populations i n Galactic star 


form ing regions, for example the MYStlX (IFeigelson et al.ll2013l) and MOXG (ITownslev et ah 


2014 ] ) projects, have led t o better censuses of young stars (both high and low-mass; both 
disk-bearing an d disk-free: iBroos et al.ll2013h that reveal snapshots of the structure of young 
stellar clusters flKnhn et al .112014 henceforth Paper I). The complex cluster structures seen 
in these regions resemble structures produced by hydrodynamical simulations of star-cluster 


form a tion from turbul e nt molecular cloud s fe.g.. lBate et al.ll2003[lBatell2009al:iBonnell et al. 


2011; 

Dale et al. 

2012; 

Walch et al. 

2012) 


formation and evolution unfold in real clusters in the Galaxy, through comparison of the 
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properties of a large sample of clusters spanning a range of environments and ages, it is pos¬ 
sible to make inferences about processes such as gas removal, cluster expansion, subcluster 
mergers, dynamical bounded ness and relax a tion, a nd cluster dispersal, which all affect the 
spatial structure of clusters. iGetman et al.l f)2014bl) provide age estimates for over 100 sub¬ 
clusters in MYStlX MSFRs, which span a wide age range between 0.5 and 5 Myr, showing 
that individual regions often have spatially segregated structures with different ages. The 
age information is used here to examine how subclusters dynamically evolve. 


The organization of this paper is as follows. Section 2 describes the data that are 
available for stars and stellar subclusters from the MYStlX project. Section 3 uses multi¬ 
variate statistics to compare measured physical properties of subclusters. Section 4 derives 
dynamical properties for the subclusters and uses multivariate statistics to investigate their 
dynamical evolution. In the conclusion to this paper (Section 5) we discuss implications for 
cluster-formation theory. 


2. Data 


This paper uses the Paper I catalog of 142 subclusters of young stars in 17 MYStlX 
star-forming regions. These subclusters were found and cha racterized us i ng the stars in the 
MYStlX Probable Complex Member (MPCM) catalogs from Broos et ah f 2013lj. The multi¬ 


wavelength data analysis ef f orts that went into thi s catalog are descri b ed bviFeigelson et ah 


( 2013). Kuhn et alJ (l2013ah iTownslev et al.l (l2014j) . King et al.l (120131) . Kuhn et al.l fl2013bl) . 


Naylor et al.l (120131 ). and IPovich et al.l (120131 ). which provided uniform data coverage across 


the 17 star-forming regions investigated here, including the most comprehensive and reliable 
lists of young stars in many of the nearest MSFRs. For the analysis of spatial structure, we 
use a subset of the MPCM sources from which weak X-ray sources have been removed in 
order to produce a sample with spatially uniform X-ray flux detection limits (Paper I). 


For cluster analysis we used the “hnite mixture model” method. To construct these 
models, the projected spatial distribution of stars in each individual subcluster is modeled 
with its own probability density, which has the form of an “isothermal ellipsoid” (Equa¬ 
tion 4 in Paper I). The isothermal ellipsoid is roughly constant within a core radius, rc, 
but decreases with a -2 power law-index outside this radius. Each subcluster is described 
by six parameters: cluster center (R.A., Dec.), subcluster core radiu^ (rc), ellipticity (e), 
orientation (0), and central surface density (Eq). In this case, the hnite mixture model is 
surface-density distribution obtained through summation of every isothermal ellipsoid for 


^The radius is the harmonic mean of the semi-major and semi-minor axes of the core ellipse. 
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all subclusters plus a model component for unclustered stars (Equation 5 in Paper I), and 
the likelihood of a particular hnite mixture model given a set of points can be calculated 
using Equation 6 in Paper I. The best-htting parameters for each subcluster are found us¬ 
ing maximum likelihood estimation, and the number of subclusters is determined using the 
Akeike Information Criterion, a likelihood that is penalized by the total number of parame¬ 
ters in the model. We find good agreement between stellar surface density estimate from this 
parametric model and surface densities estimated using an adaptive smoothing algorithm. 


Kuhn et ahl (120151 henceforth Paper II) provides estimates of the total number of stars 
in each subcluster. The number of low-mass s tars missing from the MPCM catalogs is 
inferred from the X-ray luminosity function (e.g., iGetman et al.ll2005l: iFeigelson et ahl 120051) 
and initial mass function, which provide independent but consistent estimates of intrinsic 
populations]^ Since Paper I does not provide outer truncation radii for subclusters, we 
must chose a characteristic radius to describe the size of clusters. Here, we define r 4 = 4rc 
as a characteristic radius—this coincides with our radius choice for determining subcluster 
assignments for individual stars in Paper I and would be roughly the projected half-mass 
radius for a subcluster with an outer truncation radius at ~ 17rc. We define n 4 to be the 
number of stars within a projected radius r 4 ; and, following Equation 3 of Paper I, 714 , r 4 , 
and Eq are related to each other by the equation 


774 = 0.56 rl Eq. (1) 

The conversion between the central volume density and central surface density of an isother¬ 
mal sphere is 

Po = Eo/2rc = 2Eo/r4. (2) 

And, we use the above equation as the definition of the parameter po tabulated in this 
paper—this relation between central volume density and projected surface density would be 
approximately true for the isothermal ellipsoid if the subcluster is neither much more com¬ 
pressed nor extended along the line of sight than in the plane of the skyj^ As a consequence 
of these dehnitions, any two of the variables r 4 , 774 , Eq, and po that we use to describe the 
modeled ellipsoids are sufficient to fully describe the model. 


^We choose to use the estimates obtained via the X-ray luminosity function because there are fewer 
sources of systematic uncertainty in estimated X-ray luminosities compared to estimated stellar masses used 
by Paper II. 


^We cannot avoid this assumption because we do not have measurements of the small differences in 
distance to different c luster member s. Similar approximations have been made by other studies of cluster 
central densities (e.g.. IPfalznerll2009l l. 
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The obscuration of subclusters is measured using the near-IR J — H color index and the 
X-ray median energy (ME) indicator. Paper I provides median J — H and Mi? values for the 
MPCMs assigned to each subcluster (larger values indicate higher obscuration). Although, 
obscuration is a projection effect, subclusters with greater obscuration tend to be more 
deeply embedded in the molecular cloud. 


Subcluster ages are obtained bv iGetman et al.l fj2014bl) using the novel Agejx method, in 
which X-ray luminosities of low-mass pre-main sequence stars are used as a proxies for stellar 
masses and dereddened J-band luminosities as proxies for bolometric luminosities. Ages for 
individual young stars may have large uncertainty, due to statistical error in luminosities, 
uncertainties in dereddening, and the inherent scatter in the X-ray luminosity~mass relation. 
Nevertheless, we are not so much interested in exact ages for individual objects as we are 
in consis tent treatm e nt th at permits comparison between subclusters. Ages are calculated 
using the lSiess et ahl fjl997l) PMS stellar evolution tracks—different evolutionary models may 
systematically shift ages, but are unlikely to change the ordering of the median stellar ages 
of different subclusters. 


Paper I dehned four heuristic classes of morphological structures (i.e., the arrangements 
of subclusters) seen in MSFRs, which include “simple” structures composed of an isolated 
isothermal ellipsoid, “core-halo” structure, “clumpy” structures, and long “chains” of sub¬ 
clusters. We label the subclusters by the class of the large-scale structure that they are a 
component of. For example, the star-forming region DR 21 has a chain structure, in which a 
line of subclusters is embedded in a dense molecular cloud hlament, so all DR 21 subclusters 
(A through I) are labeled with “chain.” On the other hand, some regions show multiple 
structures; for example, the Carina region has an overall chain-like structure, but the Tr 14 
cluster (our subclusters B and C) have a core-halo formation (Paper I, their Figure 2). 
These classihcations are based on the subclusters’ placements relative to other subclusters 
in a region, rather than on the subclusters’ own intrinsic properties. 

Table [T] presents the astrophysical properties of subclusters: the morphological class 
that a subcluster belongs to, subcluster size out to four core radii the number of stars 
projected within four core radii (n 4 ), the central density of stars within a subcluster (cxo 
and Pq), subcluster ellipticity (e), median interstellar medium absorption indicators {J — H 
and ME), and median stellar age {Agejx)- Logarithmic values are used for radii, numbers 
of stars, densities, and ages due to the large dynamic range of their values. We preserve 
two signihcant hgures beyond the decimal point for the logarithmic values to preserve the 
peaks in the distributions when un-logged values are used; however, uncertainties on radius, 
number of stars, and age are nearer 0.1 dex {^25%). 

In many of the scatter plots in this paper, characteristic error bars are shown as -|- 
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symbols representing the median statistical uncertainties calculated for subcluster proper¬ 
ties. The Fisher information matrices for the maximum-likelihood cluster models provide 
estimates for uncertainty on r 4 and e. These estimates do not account for systematic un¬ 
certainties in model selection but are useful for estimating the general effect of uncertainty 
on measurements. The uncertainties reported for 714 come from the scatter in number of 
observed stars when drawing from a Poisson distribution, scaled to the inferred number of 
stars. The 774 uncertainties (listed in Paper II, their Table 2) do not take into account the 
systematic errors in estimating total stellar populations from observed numbers of stars that 
arise from both the extrapolation using the X-ray luminosity function and the systematic 
uncertainty in subcluster model. Uncertainties in r 4 and 774 are propagated to Sq and po, 
assuming that uncertainty is normally distributed. Uncertainty on the medians of J — H, 
ME, and Agejx are calculated via bootstrap resampling fSee iGetman et al.ll2014bl. their 
Table 2). The resulting typical uncertainties are 0.091 dex for r 4 , 0.25 for e, 0.070 dex for 
774 , 0.19 dex for Eq, 0.283 dex for po, 0.093 dex for Agejx, 0.057 mag for median J — H, 
and 0.073 keV for median ME. Nevertheless, we estimate scatter in regressions between 
subcluster properties empirically from the data. 
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Table 1. Intrinsic Subcluster Properties 


Subcluster 

( 1 ) 

Morph. 

Glass 

( 2 ) 

logr 4 

(pc) 

(3) 

log JV 4 
(stars) 

(4) 

log Eq 

(stars pc~^) 

(5) 

log PO 

(stars pc~ll) 

( 6 ) 

e 

(7) 

J - H 
(mag) 

( 8 ) 

ME 

(keV) 

(9) 

log Age 

(yt) 

( 10 ) 

orion B 

core 

-0.71 

2.17 

3.85 

4.86 

0.3 

0.87 

1.6 

6.04 

orion C 

halo 

-0.06 

3.21 

3.58 

3.94 

0.49 

1.05 

1.6 

6.18 

orion D 

clumpy 

-0.44 

1.94 

3.07 

3.81 

0.84 

1.17 

1.4 

6.43 

flame A 

simple 

-0.31 

2.74 

3.61 

4.22 

0.37 

1.79 

2.8 

5.90 

w40 A 

simple 

-0.19 

2.48 

3.11 

3.60 

0.04 

2.10 

2.5 

5.90 

rcw36 A 

halo 

-0.24 

2.73 

3.46 

4.00 

0.33 

1.63 

2.3 

5.95 

rcw36 B 

core 

- 1.10 

1.66 

4.12 

5.51 

0.85 

2.14 

2.8 


ngc2264 A 

chain 

-0.48 

1 . 21 : 

2.42: 

3.20: 

0.14 

0.83 

1.1 

5.95 

ngc2264 B 

chain 

-0.95 

1 . 12 : 

3.28: 

4.53: 

0.44 

0.60 

1.1 


ngc2264 C 

chain 

-0.94 

1.26: 

3.39: 

4.63: 

0.11 

0.73 

1.1 


ngc2264 D 

chain 

-0.32 

1.38 

2.27 

2.89 

0.11 

0.65 

1.1 

6.50 

ngc2264 E 

chain 

-0.05 

2.16 

2.52 

2.88 

0.47 

0.63 

1.1 

6.50 

ngc2264 F 

chain 

-0.70 

1.51 

3.15 

4.15 

0.54 

0.61 

1.0 


ngc2264 G 

chain 

-0.48 

1.76 

2.98 

3.77 

0.31 

2.11 

2.3 

6.18 

ngc2264 H 

chain 

-0.16 

1.51 

2.08 

2.54 

0.19 

0.65 

1.0 


ngc2264 I 

chain 

-0.53 

1.85 

3.16 

3.99 

0.55 

1.16 

1.6 

6.18 

ngc2264 J 

chain 

-0.25 

>1.94 

>2.64 

>3.19 

0.25 

1.22 

1.7 

6.20 

ngc2264 K 

chain 

-0.17 

2.25 

2.85 

3.32 

0.55 

0.75 

1.3 

6.34 

ngc2264 L 

chain 

-0.93 

>1.33 

>3.44 

>4.66 

0.13 

1.84 

3.5 


ngc2264 M 

chain 

-0.52 

1.69 

2.98 

3.79 

0.32 

0.73 

1.2 

6.08 

rosette A 

simple 

0.61 

2.28 

1.31 

1.00 

0.28 

0.79 

1.4 


rosette B 

clumpy 

- 0.20 

1 . 00 : 

1.67: 

2.17: 

0.00 

0.73 

1.4 

6.63 

rosette C 

clumpy 

- 0.20 

>1.26 

>1.85 

>2.36 

0.71 

0.77 

1.4 

6.61 

rosette D 

clumpy 

-0.52 

1.44 

2.73 

3.55 

0.72 

1.17 

1.4 


rosette E 

clumpy 

0.54 

2.76 

1.93 

1.69 

0.08 

0.76 

1.3 

6.48 

rosette F 

clumpy 

-0.27 

0.90: 

1 . 68 : 

2.25: 

0.51 

0.83 

1.4 

6.60 

rosette G 

clumpy 

- 0.20 




0.00 

0.76 

1.0 


rosette H 

clumpy 

0.13 

1.70 

1.69 

1.86 

0.82 

0.84 

1.3 


rosette I 

clumpy 

-0.17 

> 1.11 

>1.62 

>2.09 

0.56 

0.94 

1.6 


rosette J 

clumpy 

-0.32 




0.53 

2.18 

2.4 


rosette K 

clumpy 

-0.45 




0.00 

2.00 

2.6 


rosette L 

clumpy 

0.51 

2.46 

1.69 

1.48 

0.51 

1.10 

1.5 

6.43 

rosette M 

clumpy 

0.21 

2.16 

2.00 

2.10 

0.61 

2.21 

2.3 

6.28 

rosette N 

clumpy 

-0.15 

>1.50 

>1.98 

>2.44 

0.08 

1.62 

1.9 

6.11 

rosette O 

clumpy 

-0.35 

1.26: 

2 . 22 : 

2 . 88 : 

0.75 

1.76 

2.2 

6.23 

lagoon A 

clumpy 

- 0.02 

1.93 

2.23 

2.55 

0.55 

0.84 

1.4 

6.34 

lagoon B 

clumpy 

-0.62 

2.20 

3.69 

4.62 

0.28 

1.22 

1.8 

6.15 

lagoon C 

clumpy 

-0.16 

1.93 

2.51 

2.97 

0.52 

0.85 

1.4 

6.20 

lagoon D 

clumpy 

-0.52 

1.52 

2.82 

3.65 

0.41 

0.78 

1.3 

6.25 

lagoon E 

clumpy 

0.24 

2.32 

2.10 

2.16 

0.6 

0.83 

1.3 

6.28 

lagoon F 

clumpy 

0.56 

2.86 

2.00 

1.74 

0.36 

0.80 

1.3 

6.36 

lagoon G 

clumpy 

-0.58 

1.63 

3.04 

3.92 

0.4 

0.83 

1.3 

6.34 

lagoon H 

clumpy 

-0.03 

2.36 

2.67 

3.00 

0.2 

0.79 

1.3 

6.32 

lagoon I 

clumpy 

0.12 

2.53 

2.56 

2.74 

0.1 

0.83 

1.3 

6.32 

lagoon J 

clumpy 

- 0.01 

2.09 

2.36 

2.67 

0.05 

0.87 

1.3 

6.43 

lagoon K 

clumpy 

0.13 

2.39 

2.38 

2.55 

0.45 

1.03 

1.4 

6.15 

ngc2362 A 

simple 

-0.18 

1.83 

2.44 

2.93 

0.38 

0.65 

1.2 

6.50 

ngc2362 B 

simple 

0.21 

2.60 

2.43 

2.53 

0.1 

0.62 

1.09 

6.46 

dr21 A 

chain 

-0.35 

>1.32 

>2.27 

>2.92 

0.69 

1.97 

2.3 

5.78 

dr21 B 

chain 

-0.63 

>0.89 

>2.38 

>3.31 

0.34 

2.04 

2.8 


dr21 C 

chain 

-0.64 

1.56: 

3.09: 

4.03: 

0.41 

3.00 

3.5 


dr21 D 

chain 

-0.37 

2.25 

3.25 

3.92 

0.61 

2.92 

4.0 

5.84 

dr21 E 

chain 

-0.18 

2.30 

2.92 

3.40 

0.44 

2.36 

3.7 

6.00 

dr21 F 

chain 

-0.79 

> 1.11 

>2.74 

>3.83 

0.55 

2.55 

4.0 


dr21 G 

chain 

-0.32 

>1.52 

>2.40 

>3.02 

0.22 

2.51 

3.3 


dr21 H 

chain 

-0.52 

1.77: 

3.06: 

3.87: 

0.53 

2.96 

3.3 


dr21 I 

chain 

-0.37 

>1.51 

>2.51 

>3.18 

0.32 

2.66 

3.0 

6.04 

rcw38 A 

halo 

0.86 

3.68 

2.01 

1.35 

0.19 

1.36 

2.2 

6.28 

rcw38 B 

core 

-0.53 

3.51 

4.61 

5.34 

0.36 

1.06 

2.6 

6.28 

rcw38 C 

clumpy 

- 0.22 

2.21 

2.68 

3.09 

0.78 

1.62 

2.5 

6.28 

ngc6334 A 

chain 

-0.31 

2.23 

3.11 

3.72 

0.31 

1.44 

2.1 


ngc6334 B 

chain 

0.16 

2.62 

2.56 

2.70 

0.25 

1.44 

2.0 

6.36 

ngc6334 C 

chain 

-0.64 

2.07 

3.61 

4.56 

0.15 

1.10 

1.8 


ngc6334 D 

chain 

-0.51 

2.08: 

3.35: 

4.17: 

0.09 

1.88 

2.8 


ngc6334 E 

chain 

-0.03 

2.74 

3.05 

3.37 

0.24 

1.88 

3.1 













Table 1—Continued 


Subcluster 

( 1 ) 

Morph. 

Class 

( 2 ) 

logr 4 

(pc) 

(3) 

log JV 4 
(stars) 

(4) 

log Eq 

(stars pc~^) 

(5) 

log PO 

(stars pc~ll) 

( 6 ) 

e 

(7) 

J - H 
(mag) 

( 8 ) 

ME 

(keV) 

(9) 

log Age 

(yt) 

( 10 ) 

ngc6334 F 

chain 

-0.09 

2.18 

2.61 

3.00 

0.17 

1.28 

1.8 


ngc6334 G 

chain 

-0.37 

2.41 

3.4 

4.07 

0.26 

1.61 

2.5 


ngc6334 H 

chain 

- 0.22 

2.34 

3.03 

3.55 

0.26 

1.38 

1.8 

6.20 

ngc6334 I 

chain 

-0.35 

1.91: 

2 . 86 : 

3.52: 

0.32 

1.01 

1.6 


ngc6334 J 

chain 

0.05 

3.02 

3.17 

3.42 

0.65 

2.27 

3.2 

6.18 

ngc6334 K 

chain 

-0.34 

1 . 66 : 

2.60: 

3.24: 

0.31 

3.00 

3.1 


ngc6334 L 

chain 

- 0.00 

2.51 

2.77 

3.07 

0.16 

2.44 

3.2 

5.84 

ngc6334 M 

chain 

-0.15 

>1.24 

>1.79 

>2.24 

0.2 

2.03 

2.7 


ngc6334 N 

chain 

0.05 

>1.41 

>1.56 

>1.81 

0.3 

2.20 

1.6 


ngc6357 A 

simple 

- 0.01 

3.09 

3.37 

3.68 

0.22 

1.26 

1.9 

6.15 

ngc6357 B 

simple 

0.27 

2.96 

2.67 

2.70 

0.36 

1.30 

1.9 

6.15 

ngc6357 C 

clumpy 

0.04 

3.02 

3.2 

3.46 

0.19 

1.29 

1.9 

6.08 

ngc6357 D 

clumpy 

-0.72 

2.42 

4.13 

5.16 

0.38 

1.27 

2.0 

6.04 

ngc6357 E 

clumpy 

0.24 

2 . 661 “ 

2.43^ 

2.49^ 

0.83 

1.33 

1.9 

6.15 

ngc6357 F 

simple 

0.02 

3.19 

3.40 

3.68 

0.5 

1.41 

2.1 

6.18 

eagle A 

core 

-0.33 

2.24 

3.14 

3.77 

0.01 

0.88 

1.4 

6.38 

eagle B 

halo 

0.40 

3.36 

2.81 

2.71 

0.5 

0.98 

1.5 

6.32 

eagle C 

clumpy 

-0.19 

2.07 

2.71 

3.21 

0.6 

0.96 

1.4 

6.23 

eagle D 

clumpy 

0.62 

3.15 

2.17 

1.85 

0.61 

1.10 

1.6 

6.40 

eagle E 

clumpy 

-0.53 

>1.53 

>2.84 

>3.67 

0.07 

2.40 

2.6 

6.00 

eagle F 

clumpy 

0.23 

2.32 

2.12 

2.19 

0.13 

1.61 

2.3 


eagle G 

clumpy 

-0.24 

1.75: 

2.49: 

3.03: 

0.51 

2.59 

3.7 


eagle H 

clumpy 

-0.42 

>1.35 

>2.45 

>3.17 

0.26 

1.57 

1.6 


eagle I 

clumpy 

0.02 

1.91 

2.12 

2.40 

0.64 

1.64 

2.4 

5.90 

eagle J 

clumpy 

-0.04 

1.59: 

1.92: 

2.25: 

0.41 

1.51 

2.1 


eagle K 

clumpy 

-0.61 

>1.35 

>2.82 

>3.73 

0.04 

1.42 

2.0 

6.28 

eagle L 

clumpy 

-0.92 

>1.14 

>3.25 

>4.47 

0.31 

2.01 

3.0 


ml7 A 

clumpy 

-0.57 

2 . 20 : 

3.60: 

4.47: 

0.03 

1.92 

3.2 


ml7 B 

clumpy 

- 1.02 

1.97: 

4.27: 

5.59: 

0.1 

2.49 

2.9 


ml7 C 

clumpy 

-0.26 

2.62 

3.39 

3.95 

0.09 

1.70 

2.5 

6.15 

ml7 D 

clumpy 

0.05 

3.19 

3.33 

3.58 

0.07 

1.53 

2.4 

6.04 

ml7 E 

clumpy 

-0.49 

2.23 

3.47 

4.27 

0.18 

1.56 

2.5 

6.38 

ml7 F 

clumpy 

-0.74 

2.12 

3.86 

4.90 

0.03 

1.70 

4.4 


ml7 G 

clumpy 

-0.83 

2.01 

3.92 

5.05 

0.05 

1.55 

3.9 


ml7 H 

clumpy 

-0.35 

2.62 

3.58 

4.23 

0.07 

1.25 

2.2 

6.00 

ml7 I 

clumpy 

-0.17 

2.72 

3.32 

3.80 

0.02 

1.45 

2.1 

6.15 

ml7 J 

clumpy 

-1.37 




0.02 

1.56 

2.0 


ml7 K 

clumpy 

-0.16 

2.94 

3.52 

3.98 

0.14 

1.41 

2.2 

6.00 

ml7 L 

clumpy 

-0.26 

3.13 

3.91 

4.47 

0.05 

1.82 

2.8 

6.08 

ml7 M 

clumpy 

-0.15 

2.66 

3.21 

3.66 

0.21 

1.83 

2.5 


ml7 N 

clumpy 

-0.30 

2.50 

3.35 

3.94 

0.27 

1.33 

2.1 

6.20 

ml7 O 

clumpy 

-0.50 

2.48 

3.74 

4.54 

0.18 

1.42 

2.4 

5.84 

Carina A 

clumpy 

-0.29 

2.21 

3.04 

3.63 

0.55 

1.03 

1.6 

6.45 

Carina B 

halo 

0.35 

3.39 

2.94 

2.89 

0.35 

0.97 

1.5 

6.43 

Carina C 

core 

- 0.10 

3.13 

3.58 

3.98 

0.17 

0.94 

1.4 

6.18 

Carina D 

chain 

0.23 

2.71 

2.50 

2.58 

0.5 

0.90 

1.4 

6.38 

Carina E 

chain 

0.05 

2.53 

2.70 

2.95 

0.05 

0.84 

1.4 

6.38 

Carina F 

chain 

0.31 

2.61 

2.25 

2.24 

0.36 

0.94 

1.5 

6.58 

Carina G 

chain 

1.24 

3.94t 

l.7l1‘ 

0.771“ 

0.81 

0.91 

1.5 

6.53 

Carina H 

chain 

0.02 

2.67 

2.90 

3.18 

0.23 

0.83 

1.3 

6.45 

Carina I 

chain 

-0.09 

2.34 

2.78 

3.17 

0.24 

0.80 

1.3 

6.68 

Carina J 

chain 

0.28 

2.90 

2.60 

2.62 

0.09 

0.96 

1.5 

6.36 

Carina K 

chain 

0.08 

2.66 

2.76 

2.99 

0.28 

0.89 

1.5 

6.56 

Carina L 

chain 

0.23 

2.87 

2.67 

2.75 

0.27 

0.91 

1.5 

6.43 

Carina M 

chain 

0.20 

2.75 

2.59 

2.69 

0.32 

0.95 

1.6 

6.40 

Carina N 

chain 

0.19 

2.77 

2.64 

2.75 

0.03 

1.51 

2.1 


Carina O 

chain 

-0.15 

2.62 

3.17 

3.62 

0.09 

1.15 

1.8 

6.04 

Carina P 

chain 

0.52 

3.11 

2.33 

2.10 

0.5 

0.87 

1.5 

6.62 

Carina Q 

chain 

0.21 

2.54 

2.37 

2.46 

0.35 

0.87 

1.4 

6.63 

Carina R 

chain 

0.47 

2.71 

2.02 

1.86 

0.52 

0.99 

1.6 

6.48 

Carina S 

chain 

0.13 

2.42 

2.42 

2.60 

0.41 

0.90 

1.5 

6.46 

Carina T 

chain 

0.29 

2.77 

2.44 

2.45 

0.14 

0.92 

1.5 

6.36 

trifid A 

clumpy 

0.52 

>2.04 

>1.26 

>1.04 

0.81 

1.53 

1.3 

6.28 

trifid B 

core 

-0.29 

2.48 

3.32 

3.91 

0.2 

0.92 

1.4 

6.28 

trifid C 

halo 

0.36 

2.94 

2.46 

2.4 

0.2 

0.87 

1.3 

6.28 






Table 1—Continued 


Subcluster 

( 1 ) 

Morph. 

Class 

( 2 ) 

logr 4 

(pc) 

(3) 

log JV 4 
(stars) 

(4) 

log Eq 

(stars pc~^) 

(5) 

log PO 

(stars pc~3( 

( 6 ) 

e 

(7) 

J - H 
(mag) 

( 8 ) 

ME 

(keV) 

(9) 

log Age 
(yt) 

( 10 ) 

trifid D 

simple 

0.29 

2.44 

2.10 

2.11 

0.68 

1.27 

1.4 

6.28 

ngcl893 A 

chain 

0.17 

2.49 

2.42 

2.55 

0.26 

0.81 

1.5 

6.54 

ngcl893 B 

chain 

0.11 

2.70 

2.74 

2.93 

0.21 

0.81 

1.5 

6.42 

ngcl893 C 

chain 

-0.07 

1.91: 

2.30: 

2.67: 

0.14 

0.85 

1.5 

6.50 

ngcl893 D 

chain 

-0.98 

1.34: 

3.55: 

4.83: 

0.09 

0.75 

1.4 

6.28 

ngcl893 E 

chain 

-0.45 

>1.73 

> 2.88 

>3.63 

0.67 

0.80 

1.4 


ngcl893 F 

chain 

- 0.20 

1.78: 

2.44: 

2.95: 

0.16 

0.81 

1.4 

6.32 

ngcl893 G 

chain 

-0.03 

2.29 

2.61 

2.95 

0.33 

0.89 

1.6 

6.18 

ngcl893 H 

chain 

0.07 

2.42 

2.54 

2.77 

0.1 

0.82 

1.4 

6.28 

ngcl893 I 

chain 

0.20 

2.77 

2.63 

2.74 

0.33 

0.83 

1.5 

6.45 

ngcl893 J 

chain 

- 0.88 

1.99 

4.01 

5.20 

0.37 

0.97 

1.7 

6.15 


Note. — Properties of individual ellipsoidal subclusters in Paper I. Column 1: Name of MYStlX subcluster. Column 2: 
Fiducial subcluster radius (= 4rc). Column 3: Number of stars within 4 core radii. Column 4: Central surface density. 
Column 5: Estimated central volumetric density. Column 6 : Ellipticity. Column 7: J — H NIR absorption index. Column 8 : 
X-ray median energy absorption indicator. Column 9: Median age. 
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3. Multivariate Analysis of Cluster Properties 


Figure 1—a multivariate “pairs plot”—graphically displays the astrophysical properties 
from Table 1. The rows and columns in the 9x9 array of plots corresponds to individual 
variables, whose labels and units are shown in the panels on the diagonal. The panels on the 
diagonal shows single variable distributions as histograms, the lower triangle shows bivariate 
scatterplots with LOWESS non-parametric regression lines, and the upper triangle gives 
the statistical significance of the correlations using the Kendall’s r rank test. Graphs and 
computations were performed in the R statistical software environment (R Core Team 2014). 


The LOWESS lin es show the locally-weighted polynomial regression flBecker et al.lll988 


Clevelandlll979L 119811) for the variables of the rows (dependent variable) against the variables 


of the columns (independent variable). These lines are intended to guide the eye, rather than 
be used for statistical inference, and they may be inaccurate near the edges of the distribution 
where there are fewer points. 


To test correlation between variables, our null hypothesis, Hq, is that the variables are 
uncorrelated, while our alternate hypothesis, Ha, is that a correlation exists, and p{Ho) is 
the probability of obtaining a value of r greater than or equal to the observed value. Null- 
hypothesis probabilities labeled p > 0.05 do not pass the traditional threshold for rejection 
of Hq, probabilities with 0.05 > p > 0.0001 are marginall y statist i cally s ignificant, while 
probabilities p < 0.0001 indicate strong correlations. The iKendalll (119381) rank test helps 
to reduce the effects of outlying points and uncertainties in measurement, but could still 
be vulnerable to correlated uncertainties. For example, the cluster radii and the central 
surface density (obtained from the models in Paper I) are used to calculate the numbers 
of stars and the central volume density. On the other hand, some quantities like ME and 
J — H are obtained from independent astronomical measurements, so correlations between 
these quantities and other parameters are less likely to be a result of systematic biases. For 
the hypothesis test, values are rounded to one significant figure beyond the decimal point, 
representing the approximate precision of the data. The rounding does not affect the p-values 
in an appreciable way (i.e., few ties are produced in for the rank test). 

Below, we summarize the statistical correlations with various variables that can be seen 
from Figure 1. 


Radius There is a highly statistically significant positive correlation between subcluster ra¬ 
dius and age, which could be explained astrophysically by subcluster expansion. There 
are also strong correlations with number of stars (positive), central surface density 
(negative), and central volume density (negative). Less obvious to the eye, but still 
statistically signihcant, are the negative correlations with the J — H and ME absorp- 
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tion indicators. 


Ellipticity The ellipticities measured for subclusters show little correlation with any of the 
other subcluster properties. There is a marginal indication {p < 0.05) that the denser 
subcluster are also more spherical, but there is much scatter in these relations. 

Number of Stars There is a lot of scatter in many of the correlations between the number 
of stars in a subcluster and other subcluster properties, which is somewhat surpris¬ 
ing given that this is one of the most fundamental subcluster properties. There is a 
marginally signihcant negative correlation between the number of stars and the cen¬ 
tral volume density of the cluster, which will be investigated in more depth in a later 
section. 


Density There is a tight correlation between the surface density and volume density, which 
is not a surprising result, but would not necessarily be true for every possible config¬ 
uration of subcluster structure. There is also a statistically significant trend that the 
denser subclusters have higher absorption indicated hj J — H and ME. 


J — H and ME There is a tight correlation between these two indicators of absorption by 
the molecular cloud, J — El tracing NIR dust absorption and ME tracing X-ray gas 
absorption, with a few outlying p oints. The strong inver se relation between subcluster 
age and absorption was noted bv iGetman et al.l (l2014bl l. 


Age Subcluster age has statistically significant relations with subcluster radius (positive), 
surface density (negative), volume density (negative), J — H (negative), and ME (neg¬ 
ative); which indicates that subcluster age may be an important explanatory variable 
for subcluster properties. The age of a subcluster appears to have little relation to the 
number of stars in the subcluster. 


Morphological class These classes are categorical labels without an intrinsic ordering 
(1 = “simple,” 2= “core,” 3= “halo,” 4= “clumpy,” 5= “chain”), so we do not compute 
p-values using Kendall’s r. Regarding radius, “simple” subclusters are similar in size 
to the “halo” component of core-halo structures, while the “core” components are a 
factor of ~10 smaller. The “clumpy” and “chain” subclusters cover this whole range in 
size. This trend is again reflected in the central surface and volume densities. Regard¬ 
ing the absorption indicators, J — H and ME, “simple” and “core-halo” subclusters 
typically have lower absorptions, while “clumpy” and “chain” subclusters have a wider 
range of absorptions. Regarding age, there is no obvious difference in typical age for the 
different groups; however, the “clumpy” and “chain” subclusters have a wider range of 
ages. 
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3.1. Principal Component Analysis 

Principal component analysis (PCA) is performed using the variables logr 4 , logn 4 , 
e, ME and log Age (we exclued logSo, logpo, and J — H because these quantities are 
redundant). Subclusters with missing data are excluded from this analysis. The distribution 
of ME deviates strongly from normality; most subclusters are lightly obscured. ME is 
peaked at a low value of 1.4 keV, but there is a tail of highly obscured subclusters out to 
3.0 keV. Instead we use a standardized variable based on the rank of ME values, which is 
normally distributed with a standard deviation of 0.5 (similar to other variables) as provided 
by the R-code below. 


ME.standardized <- qnorm((rank(ME)-0.5)/length(ME))/0.5 


Table [2] provides PCA loadings for the first three components. The first two principal 
components, Comp.l and Comp.2, account collectively for 84% of the variance, shown in 
the bar-chart in Figure [ 2 ] (left), so the remaining components are not particularly important. 
The components Comp.l and Comp.2 are both linear combinations of radius, number of 
stars, median energy and age, while the third principal component, Comp.3, also includes 
ellipticityfl This demonstrates that the global distribution of subcluster properties can be 
reduced down to two variables, the first one being a combination of age, absorption, and 
central (surface) density, and the second one being the number of stars in a cluster. 

Figure [2] (right) is the biplot for the first two components, Comp.l and Comp.2)] i.e., 
it shows the clu ster paramete r s proi ected into the {Comp.l, Comp.2) plane using the biplot 
definition from iGabriel et al.l (119711 1. Red arrows show the original variables projected in 
this coordinate system. The arrows for log Age, ME, and logSg are almost parallel (logpo 
points more-or-less in the same direction), which indicates that subcluster age, density, and 
absorption are all closely related. In contrast, the log 77,4 arrow is nearly orthogonal, indicating 
that—globally—the number of stars in a cluster is not strongly affected by subcluster age, 
surface density, or absorption. The r 4 arrow lies in between these two axes, reflecting our 
earlier hnding of a statistically signihcant relations between age and radius (positive) and 
radius and number of stars (positive). 


^When principal component analysis is performed using the original ME rather than the standardized 
ME, the relative relations between the variables remain the same albeit rotated in the {Comp.l, Comp.2) 
plane. 
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Table 2. PC A of Sub cluster Properties 


Property 

Variance 

Comp.l 

56% 

Comp.2 
28% 

Comp.3 
9% 

log r 4 

-0.361 

-0.471 

-0.630 

log n 4 

-0.892 


0.350 

e 



-0.549 

ME^t. 

-0.270 

0.812 

-0.384 

log Age 


-0.336 

-0.179 


-4 -2 0 2 4 6 




Component 1 

Fig. 2.— Left: A scree plot showing the relative contributions of the hve PCA components 
to the total variance. Right: A biplot shows the subcluster observations (black points) and 
the Table 1 variables (red arrows) projected onto the plane of the hrst two PCA components. 
The variables for age, absorption, and stellar density are roughly parallel, while the number 
of stars in a cluster is orthogonal, and cluster radius is in between. 



























3.2. The Subcluster Age—Radius Relation 


Figure [3] shows the relation between logr 4 and log Age. Although there is signihcant 
scatter in this relation, there is a clear, positive trend between subcluster median age from 
the AgeJX analysis and subcluster size. This can be interpreted as cluster expansio n—a phe- 
to occur as gas is removed from young embedded clusters fe.g. iTutukov 


1978; 

Hills 

1989: G 

0 

0 

Bastian 

2011 ) 


Goodwin fc Bastiaru 120061: IBaumgardt fc Kroupal 120071: iMoeckel fc Bate 


Scatter in the relation is caused by uncertainty in estimating stel¬ 
lar ages and subcluster radii, but it can also be intrinsic; i.e. subclusters form with different 
initial radii, or expansion happens at different rates for different subclusters. We £t a sim¬ 
plistic linear regression model, which helps provide intuition about this relation, even if it 
is not physically realistic. For this regression analysis, both variables have uncertainty and 
it is unclear if one variable can be labeled the “independent” variable and the other the 
“dependent ” varia ble: thus, we show the results of several linear regression strategies from 


Isobe et al.l fllOOOl) . 


The gray, dashed lines on the plot indicate the size of a sphere expanding from the center 
of the cluster at a constant velocity of 0.5, 1, 2, and 4 km s“^. Most of the regression lines 
are slightly steeper than these lines (the ordinary least squares regression does not account 
for uncertainty in median age, so its slope is likely to not be steep enough). Few points lie 
above the 1 km s“^ line on the left side of the plot, but 10 lie in this region on the right side 
of the plot. However, the OLS bisector line and the reduced major-axis line, which look like 
reasonable hts, have slopes, /3, close to the (3 = 1 constant velocity slope. The parameters 
for the reduced-major axis regression line are provided in Table [3l 


A simple model for subcluster expansion is the case where the kinetic energy of the stars 
is much greater than the potential energy of the clusters, so the stars coast outward at a 
roughly constant velocity. Many st udies of stellar velocities in young stellar clusters, both 
obser vational (Orion and NGC2264 : Fiiresz et ah 2006 . 2008 : Tobin et al.ll2015 j or theoretical 
(e.g., IBate et al.l l2003l: lBatell2009al) . show initial stellar velocities of >3 km s“h However, 
some young stellar clusters hav e lower measured vel ocity dispersions, like 7 Velorum with 0.34 
and 1.60 km s“^ components fj Jeffries et al.N2014f) . Almost all of the MYStlX subclusters 
on the plot lie below the 1 km s“^ line. If a cluster were unbound with stars moving at 
velocities >3 km s“^, it is unlikely that so many stars (the radius contains on the order 
of half the stars in the stars in the cluster) would be so near the location at which they were 
formed. This suggests that the stars are not just freely drifting, but it is likely that most of 
the subclusters are currently gravitationally boundE Gonhrmation of the bound or unbound 


®Alternatively, some of the subclusters that have small r 4 for their Agejx values, could have recently 
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state of these clusters would require radial velocity observations. 


The expansion rate of boun d young stell ar clusters was also investigated bv iMoeckel &: Bate 


fl2010l) . who take the cluster of iBatel fl2009al) and evolve it forward in time to 10 Myr under 
a variety of assumptions about gas removal timescales. They hnd (e.g., their Figure 5) that 
although the initial velocity dispersion is several km s“^, the velocity dispersion decreases 
as the cluster expands to around 1 km s“^ after 2 Myr. This is reflected in the size of the 
clusters with time (e.g., their Figures 4 and 6) where the half-mass radius typically expands 
by a factor of 10 between 1 and 10 Myr (typically at a rate of ~1 km s“^). These expansion 
rates are consistent with our observations. 


The age-radius scatter plot shows little evidence for larger radii with age until 2 Myr 
(log Age = 6.3), which is also consistent with the simulation results. The LOWESS curve 
for this relation (Figured]) also shows a hint of concavity at ~2 Myr. 


become unbound. Nevertheless, the unbinding event (by gas expulsion) would have had to happen very 
recently compared to the subclusters’ total ages for this trend to still be seen. It is unlikely that the 
majority of subclusters happen to be seen in a state immediately after gas expulsion. 
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Table 3. Relations of of Subcluster Properties 


Relation 


Parameters 

Stat. 



log A 

a 

Signif. 

(1) 


(2) 

(3) 

(4) 

log r 4 = log A + a 

log Age 

-18.0 ±0.1 

2.9 ±0.4 

<0.0001 

log So = log yf -1- a 

' log r 4 

2.6 ±0.1 

-1.9 ±0.2 

<0.0001 

log po = log A +a 

log r 4 

2.9 ±0.2 

-2.6 ±0.1 

<0.0001 

log n 4 = log A + a 

log r 4 

2.5 ±0.1 

1.7 ±0.3 

<0.0001 


Note. — Reduced-major axis regression lines. Column 1: For¬ 
mula for the relation between the two quantities. Column 2; Scale. 
Column 3: Power-law index and uncertainty. Column 4; The null- 
hypothesis probability, p{Hq), that a Kendall’s r statistic greater 
or equal to the calculated value would be produced. 
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Fig. 3.— Scatter plot of subcluster median age vs. radius, r 4 . Linear regression lines are 
shown as black lines: OLS(|/|a;) (dot-dashed), OLS(a;||/) (dashed), OLS bisector (dotted), 
orthogonal regression (solid), reduced major-axis (long dashed). Gray dashed lines indicate 
the age-radius relation for a sphere expanding from the origin at a constant velocity of 0.5, 
1, 2, and 4 km s“^. 
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3.3. The Radius—Central Density Relation 


If clusters are expanding, the density of stars in the center of the cluster will become 
attenuated with increasing radius. Figure 0] shows the surface density, Sq, radius relation 
and the volume density, po) radius relation for subclusters. Both these graphs contain the 
same information, since volume density is inferred from surface density and radius. However, 
volume density is more astrophysically meaningful, while surface density can be estimated 
without assuming a subcluster radius. The distribution of points suggests that the relation 
between radius and density is nearly a power law. The same linear regressions strategies 
that were used above are used here (the black lines on the plot). 


If a subcluster were expanding isomorphically, with no stars gained or lost, the central 
surface density would decrease proportionally to and the central volume density would 
decrease proportional to r~^ (these lines are shown in gray). However , the o bserved relations 
are somewhat less steep, with Sq oc and Pn oc iPfalznerl (1201 ll) notes a similar 


trend for young embedded clusters from ( Lada &: Lada 20031) where po 


oc r 


- 1 . 3 ± 0.7 


These flatter slopes indicate that larger clusters contain more stars. The number of 
stars within a projected radius is a function of radius and surface density, so the 714 ~ r 4 
relation shown in Figure 0 is another method of displaying the information from Figure 01 
The parameters for the reduced-major axis regression lines described above are provided in 
Table 01 
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Fig. 4. — Left: Scatterplot of central surface density, Sq, vs. subcluster radius, r4. Right: 
Scatterplot of central volume density, po, vs. subcluster radius, r^. Black lines indicate 
linear regression fits (line styles have the same meaning as in Figure E]). Gray lines indicate 
density-radius tracks for an expanding cluster with a constant number of stars. 



Fig. 5. — Scatterplot of number of stars, 714, vs. radius, r4. Black lines indicate linear 
regression hts (line styles have the same meaning as in Figure E]). 
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3.4. The Central Density—Number of Stars Relation 


From Figure [T], one can see that surface density Sq is mostly independent of the number 
of stars in a cluster while there may be a slight negative relation between po and 71 , 4 . 
To investigate this further, we partition the sample into a group with 77,4 < 200 stars (54 
subclusters) and a group with 77,4 > 200 stars (63 subclusters). Figure [6] shows the po ~ '^4 
relation, and the dashed horizontal lines indicate the median po values for the two subsamples, 
revealing that the group of subclusters with more stars tend to have lower volume densities. 

We can test the statistical signihcance of this trend using the Anderson-Darling two- 
sample test—in Figure [7] (left) the cumulative distributions of Eg of the two samples (774 < 
200 stars in black; 774 > 200 stars in gray) are compared, showing very similar distributions 
with a p-value of 0.67 (not signihcant). In contrast, for the cumulative distributions of po, 
there is a visible shift, and the Anderson-Darling test gives a p-value of 0.01 (marginally 
signihcant). 

The peak surface density can be measured independently from subcluster radius—i.e. 
Paper I (their Figures 4 and 6) hnds that there is good agreement between the central surface 
densities obtained by model htting and non-parametric smoothing—so the Eg ~ 774 relation 
may be more robust than the pg ~ 774 relation. If we consider Eg, p 4 , and r 4 as functions of 
774 , then we have Eg ( 774 ) = 2r4(774)pg(774). From Figure [5] we see that that r 4 is an increasing 
function of 774 , plus scatter, which implies that pg must have a negative relation with 774 in 
order to produce the observed hat Eg ~ 774 relation. 


An inverse trend between volume density and number of stars may seem counterintuitive. 
Nevertheless, there are well-known examples that do exhibit the trend that we observe: in 
the Orion region, the Orion Ho t Core surrounding the BN/KL object is denser than the 
much richer Trapezium cluster (jRivilla et ah 2013 ). Note that the BN/KL object is likely 
behind the ONC cluster, and only coincidentally superimposed due to the direction that we 
are observing from. 


If the pg ~ 774 relation is valid, it could be explained by the merging of subclusters. 
Consider, for example, two subclusters with the same initial radii rinit and masses M. If 
they coalesce, but without time to dynamically relax, the initial potential energy will be 
approximately the same as the hnal potential energy due to the conservation of energy (i.e., 
the kinetic energies of the stars will only be slightly perturbed if the merger happens quickly 
enough). Thus, the initial potential energy when the clusters are far away from each other 
will be 

PE = kGMyr,^, + = 2kGMyn^^, (3) 

where G is the gravitational constant and fc is a constant of order unity. The hnal potential 
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energy will be 

PE = kG{2Mf/r^^^,. 

Thus, rgnai = 2 rinit, and the ratios of hnal to initial density will be 

Pfinal/Pinit ~ (2M/M)(rinit/rfinal)^ = 1/4. 


(4) 


(5) 


A similar effect is observed by ISmith et al.l fj201l[l in simulations of interacting and merging 
stellar clusters. The encounters between subclusters create low-density halos of stars around 
the clusters which help facilitate the cluster mergers. 
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Fig. 6.— Scatterplot of central volume density vs. the number of stars in a subcluster. The 
dashed lines indicate the median po for subclusters with <200 stars and subclusters with 
>200 stars. 






Fig. 7.— Left: Cumulative distributions of Sq for subclusters with <200 stars (black) and 
>200 stars (gray). The p-value for the two-sample Anderson-Darling test with the null 
hypothesis that the distributions are the same is p = 0.67. Right: Cumulative distribution 
of po for subclusters with <200 stars (black) and >200 stars (gray). The p-value for the 
two-sample Anderson-Darling test is p = 0.01. 
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3.5. The Age—Number of Stars Relation 


When considering only the properties and age for all subclusters, as shown in the 
scatter plot in Figure | 8 ], there is no sign of correlation between subcluster age and the 
number of stars in a subcluster. In addition, the principal component analysis in Figure [2] 
shows that age and number of stars are almost orthogonal quantities. Therefore, almost any 
combination of age and number of stars in a subcluster are possible. This is despite of the 
strong age ~ relation and the strong 77.4 ~ r 4 relation. 


This lack of correlation may appear to contradict the hierarchal mergers of subclusters 
scenario that was indicated by the subcluster ellipticity distribution (Paper I), the po ~ ''^4 
relation, and the 774 ~ r 4 relation. Nevertheless, this result could still be consistent with a 
scenario in which individual subclusters grow in 77,4 over time but collections of subclusters 
from different star-forming regions with different environments do not show a trend between 
age a nd 774 . For exa mple, gas expulsion in some regions may slow their clusters’ mergers 
(e.g., lKruiissenll2012l) . so different s tar forming regions at the same age could be at different 
points in their subclusters’ growth. iFellhauer et al.l (120091) suggests that subcluster mergers 
happen quickly before gas expulsion but are impeded after gas expulsion, so subclusters may 
behave differently at the same age depending on when their gas was expelled. A possible 
example of these phenomena in MYStlX can be seen in NGC 1893—this is one of the oldest 
star-forming regions, but it is still highly substructured, possibly due to early gas expulsion 
which has inhibited dynamical evolution of the stellar population. If we examine where 
subclusters from different MSFRs lie on the diagram, we hnd that they do not share the 
same locus. For example, the NGC 6357 subclusters all lie on the upper left side of the 
diagram, while the Rosette subclusters lie on the lower right. Thus, we cannot rule out 
subcluster mergers. 


When we attempt to control for this effect by considering only sub clusters that are 
highly embedded in molecular clouds (i.e., ME > 2.5 keV) there is a marginally significant 
positive 774 ~ age relation (p < 0.05), although this is not seen in the sample of subclusters 
that are not embedded. This hints that growth of subclusters in time does occur—for the 
embedded subcluster population—even if the effect does not appear as a global trend when 
comparing different MSFRs. 
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Fig. 8.— Scatterplot of number of stars (71,4) vs. subcluster age. The gray dashed lines show 
tracks of constant population for subcluster’s containing 100%, 10%, and 1% of number of 
stars in the ONC (~3000 stars). 
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3.6. The Effect of Morphological Classes on Relations of Physical Properties 


We consider the possibility star-forming regions with a different morphology of cluster 
structure (the “morphological class” property in Tabled]) could be undergoing different types 
of dynamical evolution. Differences in young stellar cluster evolution have been reported 
in the literat ure, for e xamp le th e three distinct youn g-stellar-cluster evolutionary tracks 
identihed by iPfalznerl (120111) and iPfalzner et ahl (120121) for embedded, starburst and non- 
starburst clusters. However, our ability to investigate the effect of morphological class on 
other subcluster properties is limited by our small samples of “simple” and “core-halo” 
clusters. A single region with a “chain” or “clumpy” structure may have tens of subclusters, 
so their subclusters make up the majority objects of the catalog of 142 subclusters. On the 
other hand, as a result of small sample sizes, statistical tests using the 9 “simple” subclusters, 
7 “core” subclusters, or 6 “halo” subclusters are likely to be inconclusive. 

In Figure ini radius-density scatterplots are shown for four groups of subclusters, strat- 
ihed by the different morphological classes. All regions reveal similar negative correlations 
between these two quantities, which are all statistically signihcant. However, there is some 
difference between the slopes of power-law regressions to these data, which are provided in 
Table d] conhdence intervals (65%) on the power-law indices are calculated from 1000 boot¬ 
strap iterations. The most prominent difference between the correlations is that the “chain” 
subclusters produce a less steep slope , more similar to wh at was found by IPfalznerl (1201111 


for young embedded subclusters in the lhada fc Ladal (120031) sample. This indicates that sub¬ 


clusters in these star-forming regions, including NGC 2264, DR 21, NGC 6334, NGG 1893, 
and parts of Garina, are gaining more stars as they expand. In contrast, the subclusters in 
“clumpy” regions have a steeper relation which is consistent with expansion without gaining 
or losing stars. There are too few subclusters from the “simple” and “core-halo” categories 
to dehnitively determine which of the relations they more closely resemble—however the 
slopes of the regression-line hts more closely resemble the “clumpy” case. 


The data points in Figure [9] are marked by cluster age, using smaller circles for younger 
subclusters and larger circles for older subclusters; subclusters without reliable age mea¬ 
surements are marked with (-I-) symbols. In the panel showing “chain” subclusters, there is 
a clear age progression, with younger subclusters in the upper left of the panel, and older 
subclusters in the lower right—which is what would be expected from a simple scenario of 
subcluster expansion. This trend is less clear in the other cases, and there are hints that it 
may even be reversed for the “clumpy” subclusters. 

Figure Ho] shows the density-number-of-stars scatterplots for the different morphological 
classes. The negative trend is clearly visible for the “chain” subclusters at high statistical 
signihcance (Table Sj). The power-law index for this trend is a ~ —2, which is close to 
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the value obtained in Equation |5] from our hypothetical merger scenario. In the “simple” 
and “core-halo” cases, subclusters lie, more-or-less, within the same locus in this parameter 
space as the “chain” subclusters; however, there are too few data points to identify any trend. 
The subclusters in the “clumpy” sample are more spread out in this parameter space, and 
there is even the suggestion of the opposite trend in the opposite direction (not statistically 
significant). 

The subclusters on the Figure-ITUI plots are also marked by age, as above. On the panel 
showing “chain” subclusters, the younger stars are located in the upper left and the older 
subclusters are located in the in the lower right, which is, again, the expected results from 
the picture of subcluster expansion and hierarchical mergers that we have been developing 
in this paper. This age trend is not clearly apparent in the other panels. 

Paper I suggests that the “clumpy” subclusters may be a different type of structure 
than the other subclusters. These subclusters correspond to modes in stellar surface-density 
maps that were judged to be real by our AIC analysis, but they are, nevertheless, often part 
of large-scale clusters of stars, as is the case for many of the “clumpy” subclusters in M 17. 
Thus, it may be incorrect to directly compare their properties to properties of “simple,” 
“core-halo,” or “chain” subclusters, which are usually discrete structures well separated 
from each other. This difference may affect how the “clumpy” class of subclusters behave. 


Table 4. 
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Relations of Morphologically Stratified Subcluster Samples 


Morph. 

Parameters 

Stat. 

Class 

logT 

a 

Signif. 

(1) 

(2) 

(3) 

(4) 

Density 

vs. Radius: 

logpo = logR± 

a log r 4 

simple 

3.2 ±1.1 

-3.7: 

0.01 

core 

3.4 ±0.6 

-2.2 ±0.7 

0.05 

halo 

3.6 ±0.9 

-2.6: 

0.02 

clumpy 

2.7 ±0.5 

-3.1 ±0.3 

<0.0001 

chain 

2.9 ±0.2 

-2.0 ±0.1 

<0.0001 

Density 

vs. Number; log po = log A-\- a log n 4 

chain 

7.3 ±0.4 

-1.8 ±0.4 

<0.0001 


Note. — Reduced major-axis regression lines 
for subcluster samples stratified by the four mor¬ 
phological classes. Column 1: Morphological 
class. Column 2: Scale. Column 3: Power- 
law index and uncertainty. Column 4: The null- 
hypothesis probability, p{Hq), from Kendall’s r 
test. 
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log r4 [pc] 


Fig. 9.— Scatter plots of po vs. r 4 stratified by the four classes—from left to right and top 
to bottom: “simple,” “core-halo” (black indicates “core”; gray indicates “halo”), “clumpy,” 
and “chain.” Symbol size indicates subcluster age, with smaller circles indicating younger 
subclusters, larger circles indicating older subclusters, and +’s marking subclusters with 
missing ages. The dashed lines show the reduced major axis regression line for statistically 
signihcant relations. 
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Fig. 10.— Scatter plots of po vs. n 4 stratified by the four classes—from left to right and top 
to bottom; “simple,” “core-halo” (black indicates “core”; gray indicates “halo”), “clumpy,” 
and “chain.” Symbol size indicates subcluster age, with smaller circles indicating younger 
subclusters, larger circles indicating older subclusters, and -f-’s marking subclusters with 
missing ages. The dashed lines show the reduced major axis regression line for statistically 
signihcant relations. 
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3.7. Temporal Evolution and Timescales 


We investigate temporal evolution of the principal component of the structural prop¬ 
erties (logn 4 , logr 4 , log So, and logpo) for each morphological class. Using the relations in 
the dehnition of these quantities from Equations [1] and |2l the principal component, a linear 
combination of any two of these variables, can be expressed as 


logK = Cl logpo + C 2 logr 4 const. = {Ci - C' 2 / 3 ) logpo + (<^ 2 / 3 ) log 77,4 const. (6) 

where Ci and C 2 are the PCA loadings, which are related to the power-law indices, a, of the 
Po ~ ?"4 relations from Table 0] by 


Cl 


-|- 1 


and Co = 




+ 1 


(7) 


The large amounts of scatter in relations between subcluster properties and subcluster 
ages make it difficult to determine the best functional form when a statistically signihcant 


tial time evolution for gas expulsion (e. 


Ybarra et a 

I 2 OI 3 ) and scale-free time evolution 

for cluster expansion fe.g.. Gieles et ah 

2012: 

Pfalznei 

2011) have been explored. We also 

investigate polynomial regressions. For the regression analysis, we treat age as the indepen- 


dent variable and K as the dependent variable—the functional forms for each regression are 
shown in Equations [SHHl and the regression parameters. A, a, iV, and r are found using 
ordinary linear least squares. 


log K = log A + a log age. 

(power-law) 

(8) 

logK = log Y -I- 0.434 age/r. 

(exponential) 

(9) 


Figure im shows the various iP-variables plotted against age for each sample of subclus¬ 
ters (all subcluster, and stratihed by morphological class)—the power-law and exponential 
regression lines are superimposed. Table 0 provides the parameters for these hts, and in¬ 
formation about the correlations. In all cases a logiP ~ age has a positive relation, which 
is statistically signihcant (using Kendall’s r test) for samples including of all, “clumpy,” 
and “chain” subclusters. The coefficients of determination—a htting statistic used to 
determine what fraction of the original scatter is accounted for by the model (i.e., generally 
K, \ indicate good hts; while -C 1 indicate poor hts)—have in a similar range of 0.3- 
0.5 for both types of models, so we cannot determine whether the power-law or exponential 
regressions better described the data. However, from the K vs. age plots, there may be a 
hint of upward curvature—particularly for the “chain” subcluster sample. 
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The exponential cnrve’s e-folding timescales vary from 0.5-1 Myr, with the “clnmpy” 
snbcluster sample having the shortest r = 0.56 ± 0.14 Myr timescales, while the “chain” 
snbcluster sample having the longest r = 1.01 ± 0.2 Myr timescale. This may indicate that 
dynamical processes proceed faster in clumpy structures, where the subclusters are often 
overlapping with each other, and more slowly in chain structures, where embedded subcluster 
are often associated with a discrete molecular core. Tim escales in this range are similar to 
the 0.8 Myr e-folding timescale for gas removal found bv lYbarra et al.l (1201311 for molecular 
cl umps in the Rosette Molecu lar clouds and similar to the 1 Myr age in the numerical models 
of iGoodwin fc Bastian 
to double in size. 


200611 required for a cluster with a 10-30% star-formation efficiency 













Table 5. Fits to Structure vs. Age Relation 


Morph. 

Functional 

Parameters 


Stat. 

Class 

Form 



Signif. 

(1) 

(2) 

(3) 

(4) 

(5) 

all 

power-law 

logA = -16 ±2, a = 2.6 ±0.4 

0.32 

<0.0001 

exponential 

log A = -1.1 ±0.2, r = 0.75 ±0.11 x 10® 

0.30 

simple 

power-law 

exponential 

log A = -13:, a = 2.2 ± 1.3 
log A = -0.6:, T = 0.87 x 10®: 

0.46 

0.37 

>0.05 

core 

power-law 

log A = —21:, a = 3.4: 

0.46 

0.05 

exponential 

log A = -1.2 ±0.9, r = 0.44 x 10®: 

0.42 

halo 

power-law 

log A = —2.4:, a = 3.9: 

0.37 

>0.05 

exponential 

log A = -1.7:, r = 0.45 x 10®: 

0.29 

clumpy 

power-law 

exponential 

logA = -22 ±5, a = 3.5 ± 0.8 
log A = -1.5 ± 0.4, T = 0.56 ± 0.14 

0.36 

0.32 

<0.0001 

chain 

power-law 

exponential 

logA = -12 ±2, a = 1.8 ±0.3 
log A = -0.8 ± 0.2, r = 1.01 ± 0.20 x 10® 

0.36 

0.38 

<0.0001 


Note. — Regression fits to the relation between the principal component of the 
subcluster structure properties, K, and subcluster median age. Column 1; The sample 
for which K is derived and the regression is performed. Column 2: The equation used to 
model the K ~ age relation. Column 3: the values of the regression parameters (using 
the variable definitions in Equations [8] and [9]) and their uncertainties (assuming the 
functional form) calculated using 1000 iterations of bootstrap resampling. Column 4: 
The B? coefficient of determination. Column 5: The null-hypothesis probability, p{Hq), 
from Kendall’s r test. Values marked with a colon are highly uncertain. 
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Fig. 11.— Subcluster age vs. the principal component, K, of the subcluster structure pa¬ 
rameters, shown for the sample of all subclusters and subsamples of each morphological 
class. The y-axes show K computed for each subsample. (The po and quantities in the y- 
axis label equations are logarithmic.) Power-law (dashed lines) and exponential (solid lines) 
regressions are shown. 
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4. Dynamical Subcluster Properties 

A variety of dynamical subcluster properties can be calculated from the physical proper¬ 
ties given in Table [1] For these calculations we will ignore the effects of gas mass, considering 
only the effect of stars d own t o 0.1 Mq. The average stellar mass for stars in this mass range 
from the iMaschbergerl (120131) initial mass function is m = 0.5 Mq, which we will use to 
convert from star counts to solar mass units@ 


The free-fall time for a sub cluster is 


tff - 


Ott 


( 10 ) 


32 Gmpo 

If a subcluster is in virial equilibrium then the velocity dispersion of the stars will be 

= ^/4:7iGr‘^mpo/9. ( 11 ) 

The crossing time for a cluster requires an assumption for the veloc ity of stars; two reasonable 
possibilities include a velocity dispersion of 3 km s“^ as seen by iFiiresz et ahl (120061.120081) 
or the virial velocity dispersion from Equation [TTl These two cluster crossing times are 


Aross,i ~ r^/s km s ^ or fcrossg ~ r^/a^. 
The sub cluster relaxation time will be 

77-4 

Aelax ~ ArossTTi ; 

8 In 77-4 

and the age of the young stellar cluster in units of freiax will be 

-^relax ^-S^/Aelax- 


( 12 ) 


(13) 


(14) 


Both Equations [13] and [TT] require a choice of which formula to use for cluster crossing time. 
The dynamical quantities obtained from these equations are tabulated in Table [6] and their 
univariate distributions are shown in Figure [T2j 

The free-fall times for subclusters tend to be ~1 Myr, slightly younger than the typical 
age of the su bclusters. Thus, if star-formation happens on the one free-fall timescale as 
postulated bv lElmegreenI (120001 ). it should have already ended in most of the subclusters. 

The virial equilibrium velocity dispersions are <1 km s“^ which is significantly less than 
the global velocity dispersions seen in MSFRs like the ONC and NGC 2264. However, gas 


^Recall that our definition of po has units [stars pc ^], not [Mq pc ^]. 
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mass unaccounted for in our derivation of may raise these values somewhat. Given that 
the survival of the subclusters indicates that they are gravitationally bound, it is unlikely 
that the true intracluster velocity dispersion is too much larger than the true values when 
gas mass is accounted for. Cluster crossing times assuming ~3 km s“^ velocities are much 
less than subcluster ages, so residual structure from star formation should be dynamically 
erased, on the other hand, this will not be the case assuming velocities similar to 


The subcluster relaxation times are mostly much greater than the subcluster ages, which 
is a common result in studies of young stellar clusters in star-forming regions. Nevertheless, 
10%-40% of the subclusters have survived for several dynamical timescales. These clusters 
tend to have both small sizes and small numbers of stars, since both r 4 and U 4 increase the 
relaxation time scale in Equations [T2] and [T31 The youngest, embedded clusters often have 
low values of both and 77 , 4 , so they often have the shortest dynamical relaxation times. If 
young stellar clusters are built up by subcluster mergers, the low-n 4 clusters may dynamically 
relax before merging into larger-scale stru ctures, which might inherit the dynamically relaxed 
state of their composite subclusters fe.g.. lAllison et al.ll2009l) . 


Figure dn] shows the relation between absorption (J — H) and relaxation time, indicating 
that both lightly and heavily absorbed subclusters can have shorter relaxation times, while 
moderately embedded subclusters have longer relaxation times. 
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Table 6. Dynamical Subcluster Properties 


Subcluster 

(1) 

logtff 

(yr) 

(2) 

log Oyirial 

(km s“l) 

(3) 

log Across, 1 

(y) 

(4) 

log Across,2 

(yr) 

(5) 

log ^relax, 1 

(y) 

(6) 

log trelax,2 

(y) 

(7) 

log A^relax, 1 
(^relax,1) 

(8) 

log -lVjelax,2 
(^relax,2) 

(9) 

orion A 

4.39 

0.07 

4.17 

4.58 

4.35 

4.75 

1.80 

1.40 

orion B 

5.10 

0.01 

4.82 

5.29 

5.39 

5.86 

0.65 

0.18 

orion C 

5.56 

0.20 

5.47 

5.75 

6.91 

7.19 

-0.73 

-1.01 

orion D 

5.62 

-0.25 

5.09 

5.82 

5.48 

6.20 

0.96 

0.23 

flame A 

5.42 

0.09 

5.22 

5.61 

6.26 

6.65 

-0.35 

-0.74 

w40 A 

5.73 

-0.10 

5.34 

5.92 

6.16 

6.74 

-0.26 

-0.84 

rcw36 A 

5.53 

0.05 

5.29 

5.72 

6.32 

6.75 

-0.37 

-0.80 

rcw36 B 

4.77 

-0.05 

4.43 

4.96 

4.61 

5.14 



ngc2264 A 

5.93 

-0.59 

5.05 

6.12 

4.92 

5.98 

1.04 

-0.03 

ngc2264 B 

5.26 

-0.40 

4.58 

5.46 

4.38 

5.26 



ngc2264 C 

5.21 

-0.34 

4.59 

5.41 

4.49 

5.30 



ngc2264 D 

6.08 

-0.59 

5.21 

6.28 

5.19 

6.25 

1.32 

0.25 

ngc2264 E 

6.09 

-0.33 

5.48 

6.28 

6.04 

6.85 

0.46 

-0.34 

ngc2264 F 

5.45 

-0.33 

4.83 

5.65 

4.90 

5.71 



ngc2264 G 

5.64 

-0.31 

5.05 

5.84 

5.30 

6.09 

0.88 

0.09 

ngc2264 H 

6.26 

-0.60 

5.37 

6.45 

5.44 

6.52 



ngc2264 I 

5.53 

-0.25 

5.00 

5.73 

5.32 

6.04 

0.85 

0.13 

ngc2264 J 



5.28 






ngc2264 K 

5.87 

-0.23 

5.36 

6.06 

5.99 

6.70 

0.35 

-0.35 

ngc2264 L 



4.60 






ngc2264 M 

5.63 

-0.33 

5.01 

5.82 

5.21 

6.02 

0.87 

0.06 

rosette A 

7.03 

-0.60 

6.14 

7.22 

6.80 

7.88 



rosette B 

6.44 

-0.83 

5.33 

6.64 

5.07 

6.37 

1.57 

0.26 

rosette C 



5.33 






rosette D 

5.75 

-0.46 

5.02 

5.95 

5.03 

5.97 



rosette E 

6.68 

-0.33 

6.08 

6.88 

7.13 

7.93 

-0.65 

-1.46 

rosette F 

6.40 

-0.86 

5.27 

6.60 

4.94 

6.28 

1.66 

0.33 

rosette G 



5.33 






rosette H 

6.59 

-0.65 

5.66 

6.79 

5.87 

7.00 



rosette I 



5.36 






rosette J 



5.22 






rosette K 



5.08 






rosette L 

6.79 

-0.46 

6.04 

6.98 

6.84 

7.78 

-0.41 

-1.35 

rosette M 

6.48 

-0.46 

5.74 

6.67 

6.30 

7.23 

-0.02 

-0.96 

rosette N 



5.38 






rosette O 

6.09 

-0.63 

5.18 

6.28 

5.07 

6.18 

1.16 

0.05 

lagoon A 

6.25 

-0.46 

5.51 

6.45 

5.88 

6.82 

0.46 

-0.48 

lagoon B 

5.22 

-0.03 

4.91 

5.41 

5.50 

6.00 

0.64 

0.14 

lagoon C 

6.04 

-0.39 

5.37 

6.24 

5.75 

6.62 

0.45 

-0.42 

lagoon D 

5.70 

-0.41 

5.01 

5.90 

5.08 

5.97 

1.17 

0.28 

lagoon E 

6.45 

-0.39 

5.77 

6.64 

6.46 

7.33 

-0.18 

-1.05 

lagoon F 

6.66 

-0.28 

6.09 

6.85 

7.23 

7.99 

-0.87 

-1.63 

lagoon G 

5.56 

-0.33 

4.95 

5.76 

5.11 

5.92 

1.23 

0.43 

lagoon H 

6.03 

-0.24 

5.50 

6.22 

6.23 

6.95 

0.10 

-0.62 

lagoon I 

6.15 

-0.23 

5.65 

6.35 

6.51 

7.21 

-0.19 

-0.89 

lagoon J 

6.19 

-0.39 

5.52 

6.39 

6.03 

6.89 

0.40 

-0.46 

lagoon K 

6.25 

-0.31 

5.66 

6.45 

6.41 

7.19 

-0.27 

-1.05 

ngc2362 A 

6.06 

-0.43 

5.35 

6.26 

5.65 

6.56 

0.85 

-0.05 

ngc2362 B 

6.26 

-0.24 

5.74 

6.46 

6.66 

7.38 

-0.20 

-0.91 

dr21 A 



5.18 






dr21 B 



4.90 






dr21 C 



4.89 






dr21 D 

5.57 

-0.12 

5.16 

5.76 

5.80 

6.40 

0.05 

-0.55 

dr21 E 

5.83 

-0.19 

5.35 

6.02 

6.03 

6.70 

-0.03 

-0.70 

dr21 F 



4.74 






dr21 G 



5.22 






dr21 H 

5.59 

-0.29 

5.02 

5.79 

5.27 

6.04 



dr21 



5.16 






rcw38 A 

6.69 

-0.18 

6.39 

7.05 

8.24 

8.90 

-1.96 

-2.62 

rcw38 B 

4.70 

0.43 

5.00 

5.05 

6.70 

6.75 

-0.42 

-0.47 

rcw38 C 

5.82 

-0.38 

5.31 

6.17 

5.91 

6.77 

0.37 

-0.49 

ngc6334 A 

5.67 

-0.16 

5.22 

5.86 

5.84 

6.48 



ngc6334 B 

6.18 

-0.20 

5.69 

6.37 

6.63 

7.31 

-0.27 

-0.95 

ngc6334 C 

5.25 

-0.08 

4.89 

5.44 

5.38 

5.93 



ngc6334 D 

5.44 

-0.14 

5.02 

5.64 

5.52 

6.14 














Table 6—Continued 


Subcluster 

(1) 

logtff 

(yr) 

(2) 

log Oyirial 

(km s-1) 

(3) 

log Across, 1 

(y) 

(4) 

log tcross,2 

(yr) 

(5) 

log ^relax, 1 

(y) 

(6) 

log trelax,2 

(y) 

(7) 

log A^relax, 1 
(^relax,1) 

(8) 

log JVrelax,2 
(^relax,2) 

(9) 

ngc6334 E 

5.84 

-0.05 

5.51 

6.03 

6.54 

7.07 



ngc6334 F 

6.03 

-0.30 

5.44 

6.22 

6.02 

6.80 



ngc6334 G 

5.49 

-0.05 

5.17 

5.69 

5.93 

6.45 



ngc6334 H 

5.75 

-0.15 

5.31 

5.95 

6.02 

6.65 

0.18 

-0.45 

ngc6334 I 

5.77 

-0.31 

5.18 

5.96 

5.54 

6.33 



ngc6334 J 

5.81 

0.05 

5.58 

6.01 

6.85 

7.28 

-0.68 

-1.11 

ngc6334 K 

5.91 

-0.43 

5.19 

6.10 

5.37 

6.28 



ngc6334 L 

5.99 

-0.18 

5.53 

6.19 

6.38 

7.03 

-0.53 

-1.19 

ngc6334 M 



5.38 






ngc6334 N 



5.58 






ngc6357 A 

5.69 

0.12 

5.52 

5.88 

6.86 

7.22 

-0.71 

-1.07 

ngc6357 B 

6.18 

-0.09 

5.80 

6.37 

7.02 

7.59 

-0.88 

-1.45 

ngc6357 C 

5.79 

0.06 

5.57 

5.99 

6.84 

7.26 

-0.76 

-1.18 

ngc6357 D 

4.95 

0.14 

4.81 

5.14 

5.58 

5.92 

0.46 

0.12 

ngc6357 E 



5.78 






ngc6357 F 

5.69 

0.15 

5.55 

5.88 

6.97 

7.30 

-0.80 

-1.12 

eagle A 

5.64 

-0.15 

5.21 

5.84 

5.83 

6.46 

0.55 

-0.08 

eagle B 

6.17 

0.04 

5.94 

6.37 

7.50 

7.94 

-1.18 

-1.62 

eagle C 

5.92 

-0.30 

5.34 

6.12 

5.83 

6.61 

0.40 

-0.38 

eagle D 

6.60 

-0.17 

6.15 

6.80 

7.53 

8.18 

-1.13 

-1.78 

eagle E 



5.01 






eagle F 

6.43 

-0.39 

5.76 

6.63 

6.45 

7.32 



eagle G 

6.01 

-0.44 

5.29 

6.21 

5.54 

6.45 



eagle H 



5.11 






eagle I 

6.33 

-0.49 

5.55 

6.52 

5.92 

6.88 

-0.01 

-0.98 

eagle J 

6.40 

-0.62 

5.50 

6.60 

5.62 

6.72 



eagle K 



4.93 






eagle L 



4.61 






ml7 A 

5.29 

-0.05 

4.96 

5.49 

5.55 

6.08 



ml7 B 

4.73 

0.06 

4.51 

4.92 

4.92 

5.33 



ml7 C 

5.55 

0.00 

5.27 

5.74 

6.21 

6.68 

-0.06 

-0.54 

ml7 D 

5.74 

0.13 

5.59 

5.93 

7.00 

7.35 

-0.96 

-1.31 

ml7 E 

5.39 

-0.07 

5.04 

5.59 

5.66 

6.21 

0.72 

0.17 

ml7 F 

5.08 

-0.00 

4.79 

5.27 

5.32 

5.80 



ml7 G 

5.00 

-0.02 

4.70 

5.20 

5.14 

5.64 



ml7 H 

5.41 

0.05 

5.18 

5.61 

6.11 

6.54 

-0.11 

-0.54 

ml7 I 

5.63 

0.01 

5.36 

5.82 

6.38 

6.85 

-0.24 

-0.70 

ml7 J 



4.16 






ml7 K 

5.54 

0.12 

5.37 

5.73 

6.58 

6.94 

-0.58 

-0.94 

ml7 L 

5.29 

0.26 

5.27 

5.49 

6.64 

6.86 

-0.56 

-0.78 

ml7 M 

5.70 

-0.03 

5.38 

5.89 

6.35 

6.86 



ml7 N 

5.56 

-0.04 

5.24 

5.75 

6.07 

6.59 

0.13 

-0.38 

ml7 O 

5.26 

0.06 

5.03 

5.45 

5.86 

6.28 

-0.01 

-0.43 

Carina A 

5.71 

-0.19 

5.24 

5.90 

5.84 

6.50 

0.61 

-0.06 

Carina B 

6.08 

0.08 

5.88 

6.28 

7.48 

7.87 

-1.04 

-1.44 

Carina C 

5.54 

0.18 

5.44 

5.73 

6.81 

7.11 

-0.63 

-0.93 

Carina D 

6.24 

-0.20 

5.76 

6.43 

6.77 

7.44 

-0.39 

-1.06 

Carina E 

6.05 

-0.19 

5.58 

6.25 

6.44 

7.11 

-0.06 

-0.73 

Carina F 

6.41 

-0.29 

5.84 

6.60 

6.76 

7.53 

-0.18 

-0.95 

Carina G 



6.77 






Carina H 

5.94 

-0.11 

5.55 

6.13 

6.53 

7.11 

-0.08 

-0.67 

Carina I 

5.94 

-0.22 

5.44 

6.13 

6.14 

6.84 

0.54 

-0.16 

Carina J 

6.22 

-0.12 

5.81 

6.41 

6.99 

7.59 

-0.63 

-1.23 

Carina K 

6.03 

-0.14 

5.61 

6.23 

6.58 

7.20 

-0.02 

-0.64 

Carina L 

6.15 

-0.11 

5.76 

6.35 

6.90 

7.49 

-0.47 

-1.06 

Carina M 

6.18 

-0.16 

5.74 

6.38 

6.78 

7.42 

-0.38 

-1.02 

Carina N 

6.15 

-0.15 

5.72 

6.34 

6.78 

7.41 



Carina O 

5.71 

-0.05 

5.38 

5.91 

6.31 

6.84 

-0.27 

-0.80 

Carina P 

6.47 

-0.14 

6.05 

6.67 

7.41 

8.02 

-0.79 

-1.40 

Carina Q 

6.30 

-0.27 

5.75 

6.49 

6.62 

7.37 

0.01 

-0.73 

Carina R 

6.60 

-0.32 

6.00 

6.79 

7.01 

7.80 

-0.53 

-1.32 

Carina S 

6.23 

-0.29 

5.66 

6.42 

6.43 

7.20 

0.03 

-0.74 

Carina T 

6.30 

-0.20 

5.83 

6.50 

6.89 

7.56 

-0.53 

-1.20 

trifid A 



6.05 






trifid B 

5.57 

-0.05 

5.24 

5.77 

6.07 

6.59 

0.21 

-0.31 






Table 6—Continued 


Subcluster 

(1) 

log tjf 
(yr) 

(2) 

log (Tvij-ial 
(km s“^) 

(3) 

log tcross, 1 

(y) 

(4) 

log tcross,2 

(yr) 

(5) 

log ^relax, 1 

(y) 

(6) 

log trelax,2 

(y) 

(7) 

log Nj-elax, 1 
(^relax,1) 

(8) 

log Nj.elax,2 
(^relax,2) 

(9) 

trifid C 


6.33 

-0.15 

5.90 

6.52 

7.10 

7.73 

-0.82 

-1.45 

trifid D 


6.47 

-0.36 

5.83 

6.67 

6.61 

7.45 

-0.33 

-1.17 

ngcl893 

A 

6.25 

-0.27 

5.70 

6.45 

6.53 

7.28 

0.01 

-0.73 

ngcl893 

B 

6.06 

-0.14 

5.64 

6.26 

6.65 

7.26 

-0.23 

-0.85 

ngcl893 

C 

6.19 

-0.45 

5.46 

6.39 

5.82 

6.75 

0.68 

-0.24 

ngcl893 

D 

5.11 

-0.28 

4.55 

5.30 

4.50 

5.25 

1.78 

1.03 

ngcl893 

E 



5.08 






ngcl893 

F 

6.05 

-0.44 

5.33 

6.25 

5.59 

6.51 

0.73 

-0.19 

ngcl893 

G 

6.05 

-0.27 

5.50 

6.25 

6.16 

6.91 

0.01 

-0.74 

ngcl893 

H 

6.14 

-0.26 

5.60 

6.34 

6.37 

7.11 

-0.09 

-0.83 

ngcl893 

I 

6.16 

-0.15 

5.73 

6.35 

6.79 

7.41 

-0.34 

-0.97 

ngcl893 

J 

4.93 

0.00 

4.65 

5.12 

5.08 

5.55 

1.07 

0.60 


Note. — Dynamical properties of individual ellipsoidal subclusters. Column 1: Subcluster designation. Column 2: Free-fall time 
for the subcluster (stellar mass only). Column 3: Mean velocity dispersion for stars within 4 core radii assuming that the subcluster 
is virialized. Column 4: Subcluster crossing time assuming a velocity dispersion of 3.0 km s~^. Column 5: Subcluster crossing time 
assuming a virial velocity dispersion. Column 6: Dynamical relaxation time assuming a velocity dispersion of 3.0 km s~^. Column 7: 
Dynamical relaxation time assuming a virial velocity dispersion. Column 8—9: Subcluster dynamical age in relaxation times. 
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Fig. 13.— Scatter plot of treiax,i vs J — H, with a LOWESS line plotted. (The relation using 

treiax,2 is similar.) 
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5. Conclusions 

Our primary empirical results and their implications for star-formation theory are listed 
below. 


1. Subclusters with larger ages have larger radii, which we validate using Kendall’s r 
test (p < 10" ^) ■ We interpre t this to be a sign of subcluster expansion, which was 
postulated by iTutukovI fll978l) and others as a consequence of the loss of molecular- 
cloud material due to star-formation feedback. This helps conhrm the inference of 


cluster expansion bv lPfalzneii fcOllh from the d ensity-radius relation in the catalog of 
young stellar clusters from iLada fc Ladal (1200311 . 


The inferred expansion rate, which for most subclusters is a radial increase less than 
1 km s“^, is lower than the radial velocities observed in many spectroscopic studies of 


Fiiresz et ah 

2006, 

2008; 

Tobin et ah 

2015) 


velocity dispersions (e.g., I Jeffries et al.N2014il . If the MYStlX regions have velocity 


dispersions >3 km s“^ the subclusters would need to be gravitationally bound for 
their structure to be preserved. However, for clusters that have undergone expansion, 
~1 km s“^ expansion rates are consistent with the decrease in veloci ty dispersion 


seen m 


young stellar cluster N-body simulations jMoeckel &: Bate 201oh . This result 
has implications for the model of cluster survival by Kruiissenr (2012') because the 
fraction of stars that form in naturally bound subclusters is an intermediate step in 
determining what fraction of stars will remain in bound clusters after the end of star 
formation. Their model can be applied before the end of star formation in a region, so, 
by combining the MYStlX results with additional studies of the molecular clouds in 
these regions, the cluster-survival model can be tested on Galactic star-forming regions. 


We observe the strong correlation from iPfalznerl (120111) between cluster radius and 
cluster density. (There is a slight difference in the quantities being measured because 
they use total cluster radius while we use subcluster core radius.) We observe a similar 
trend, po ot a track slightly less steep than would be expected for expansion from 
a uniform initial state with a conserved nu mber of s tars. The flatter slope indicates 
that one of those conditions must not hold. iPfalzneii (j201lf) explains this phenomenon 
as an effect of star formation progressing from the inner part of the clusters, first, 
to the outer part o f the cluster, later, (t he inside-out scenario). However, the age 
gradients found by iGetman et ahl (j2014al) show the opposite trends, ruling out the 
inside-out scenario. We note that hierarchical cluster mergers also offer an explanation 
for this effect because, as clusters have time to expand, they have time to merge with 
























































subclusters; and mergers with subclusters can cause expansion due to conservation of 
energy. 


4. Larger subcluster central volume density, po, is correlated with smaller numbers of 
stars, 77-4 {p < 0.05; marginally signihcant)—which is consistent with our observation 
of a flat So ~ 774 relation. This effect, combined with the distrib ution of moderately 
elliptical subclusters from Paper I (cf. iMaschberger et al.l |2010[) and the less steep 
Po ~ relati on, suggests that youn g stellar clusters grow from hierarchical subcluster 
mergers fe.g.. McMillan et al.ll2007l) . If subcluster mergers are an important aspect of 
early youn g stellar cluster dyna.mical evolution, it wo uld lead to more rapid dynamical 
relaxation flAllison et al.ll2009t iParker fc Meverll2012l ). which cou ld be used to exp lain 
unexpected results, like the mass segregation seen in W 40 (eg., iKuhn et al.ll2010l) or 
well-ht isothermal sphere prohles from Paper I. 


5. The morphological classihcations of structure of star-forming regions from Paper I 
(“simple,” “core-halo,” “clumpy,” “chain”) are not correlated with age. This may 
be considered evidence against an evolutionary progression. However, an alternate 
explanation for the lack of a trend could be differences in rates of dynamical evolution in 
different enviro nments. Gas expulsion is expected to slow mergers of subclusters (e.g.. 


Kruiissed 120121) . so when “chain” clusters —which are h ighly substructured suggesting 
an early stage of evolutionary progression fjParkerl 120141) — lose their gas, as in the case 
of NGC 1893, they may maintain their initial structure inherited from the molecular 
cloud as they age. 


6 . When comparing the subcluster age and the number of stars in a subcluster for all 
MYStlX regions simultaneously, there is no evidence of a trend. However, this may be 
influenced by differing cloud environments, as we suggest for the morphological classes. 
When we control for absorption, selecting only the most absorbed subclusters for ~ 
age analysis, the number of stars does increase with radius (p < 0.05; marginally 
signihcant). Thus, the age-number of stars relation is also consistent with a picture of 
subcluster mergers. 


7. Two-body relaxation times calculated for subclusters tend to be much longer than 
the age of the subclusters. However, a fraction of the subclusters with low -774 and 
low-r 4 , have had time to dynamically relax. This is an import ant aspect of mod¬ 


els of accelerated dyn amical relaxation through cluster mergers flAllison et al.l 12009 


Parker &: Meverl 120121) . 


Finally, it is important to note that the analysis of correlations in cluster properties 
performed in this paper is sensitive to the way in which subclusters are dehned, which is 
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true for the Paper I sub cluster catalog as well as the iLada fc Ladal (120031) catalog and the 


Gutermuth et al.l (120091) catalog. For the type of analysis performed here, the MYStlX sub¬ 


cluster catalog has several advantages over these other catalogs. First, the MYStlX MPCM 
catalog contains relatively large samples of low-mass young stars (both d isk-bearing and 


disk- f ree) without cont aining too many non-cluster member contaminants (IFeigelson et ah 


20131: iBroos et al.ll2013l ). Second, the subcluster analysis avoids combining distinct clumps of 


stars into artihcial clusters, and instead focuses on the properties of subclusters. Third, para¬ 
metric modeling disentangles subcluster properties like density and radius, instead of using 
a surface-density threshold to dehne subcluster boundaries, which would lead to subclusters 
with higher surface density having artihcially larger radii (Paper I). 
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